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Abstract 

1 Although it is well known that nitrogen (N) additions strongly affect ecto- 
mycorrhizal (EM) fungal community composition, less is known about how 
different N application rates and duration of N additions affect the func- 
tional role EM fungi play in the forest N cycle. 

2 We measured EM sporocarp abundance and species richness as well as deter- 
mined the ^'^N in EM sporocarps and tree foliage in two Pinus sylvestris for- 
ests characterized by short- and long-term N addition histories and multiple 
N addition treatments. After 20 and 39 years of N additions, two of the 
long-term N addition treatments were terminated, thereby providing a 
unique opportunity to examine the temporal recovery of EM sporocarps 
after cessation of high N loading. 

3 In general, increasing N availability significantly reduced EM sporocarp pro- 
duction, species richness, and the amount of N retained in EM sporocarps. 
However, these general responses were strongly dependent on the application 
rate and duration of N additions. The annual addition of 20 kg-N ha^^ for 
the past 6 years resulted in a slight increase in the production and retention 
of N in EM sporocarps, whereas the addition of 100 kg-N-ha^'-yr^^ during 
the same period nearly eliminated EM sporocarps. In contrast, long-term 
additions of N at rates of ca. 35 or 70 kg-N-ha^'-yr^' for the past 40 years 
did not eliminate tree carbon allocation to EM sporocarps, although there 
was a decrease in the abundance and a shift in the dominant EM sporocarp 
taxa. Despite no immediate recovery, EM sporocarp abundance and species 
richness approached those of the control 20 years after terminating N addi- 
tions in the most heavily fertilized treatment, suggesting a recovery of carbon 
allocation to EM sporocarps after cessation of high N loading. 

4 Our results provide evidence for a tight coupling between tree carbon allo- 
cation to and N retention in EM sporocarps and moreover highlight the 
potential use of (5^^N in EM sporocarps as a relative index of EM fungal 
sink strength for N. However, nitrogen additions at high dosage rates or 
over long time periods appear to disrupt this feedback, which could have 
important ramifications on carbon and nitrogen dynamics in these forested 
ecosystems. 



Introduction 

It is well understood that mycorrhizas, the symbiotic 
associations between fungi and plant roots, are ubiquitous 
features of nearly all terrestrial plant communities (Smith 
and Read 2008). Mycorrhizal fungi efficiently explore the 
soil to acquire nutrients, some of which is transferred to 



host plants in exchange for recently fixed carbon (C). In 
forested ecosystems where nitrogen (N) is a major limit- 
ing nutrient, such as boreal forests, it is thought that 
ectomycorrhizal (EM) fungi play a key role for host plant 
N acquisition and tree growth (Read et al. 2004; Smith 
and Read 2008; but see Nasholm et al. 2013). However, 
increasing use of fertilizers and anthropogenic N 
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deposition from agricultural and industrial emissions has 
drastically affected the global N cycle (Vitousek et al. 
1997; Galloway et al. 2008), and, as a result, many ecosys- 
tems previously characterized by N limitation have expe- 
rienced an increase in N availability. Increasing N 
availability in natural ecosystems has major ecological 
implications, including, but not limited to, changing bio- 
logical diversity, acidifying soils, altering the global C 
cycle, and changing the rate and pathways of N cycling 
(Tilman 1987; Aber et al. 1995; Vitousek et al. 1997; 
Chung et al. 2007; Galloway et al. 2008). 

Trees decrease belowground C allocation in response to 
increasing soil N availability (Agren and Franklin 2003; 
Miikela et al. 2008; Franklin et al. 2009), which results in 
reduced C supply to EM fungi in particular (Wallenda 
and Kottke 1998; Hogberg et al. 2003; Treseder 2004; 
Hogberg et al. 2010). Numerous field studies have shown 
a reduction in C allocation to EM fungi in response to 
elevated N deposition and N additions (Arnolds 1991; 
Wallander and Nylund 1992; Termorshuizen 1993; Arne- 
brant 1994; Wallenda and Kottke 1998; Lilleskov et al. 
2001; Nilsson and Wallander 2003; Nilsson et al. 2005, 
2007; Ostonen et al. 2011; KjoUer et al. 2012), with 
changes in aboveground sporocarp production being 
more sensitive to N supply than changes belowground. 
It is also well recognized that increases in soil N availabil- 
ity reduce EM species richness and cause significant 
changes in EM fungal community composition (Lilleskov 
et al. 2001, 2002a,b; Avis et al. 2003; Toljander et al. 
2006; Parrent and Vilgalys 2007; Cox et al. 2010; Lilleskov 
et al. 2011; KjoUer et ak 2012). In short, the results show 
an increase in Lactarius, Laccaria, and Paxillus species 
with increasing N additions (i.e., "nitrophilic" taxa sensu 
Lilleskov et al. 2001), whereas Cortinarius, Suillus, Tricho- 
loma, and Piloderma species respond negatively to ele- 
vated N (i.e., "nitrophobic" taxa sensu Lilleskov et al. 
2001). Despite our increasing knowledge of the impacts N 
additions have on EM fungal community composition, 
less is known about how different N application rates and 
duration of N additions affect the functional role EM 
fungi play in the forest N cycle. Furthermore, the recovery 
of EM sporocarps after the termination of long-term N 
additions remains unclear. 

Because of the tight coupling between C and N in 
mycorrhizal symbioses, changes in tree C allocation to 
EM fungi can have dramatic effects on the transfer of N 
in EM symbioses (Alberton et al. 2007; Niisholm et al. 
2013). Although the common dogma is that N is readily 
transferred from EM fungi to host plants in environments 
characterized by low N availability, several pot studies 
(Colpaert et al. 1996; Alberton et al. 2007; Correa et al. 
2008; Alberton and Kuyper 2009) as well as a recent field 
study (Nasholm et al. 2013) have demonstrated that EM 



fungi retain large amounts of N. Using a N tracer tech- 
nique, Nasholm et al. (2013) showed that greater C allo- 
cation to EM fungi in an N-limited boreal forest resulted 
in increased N retention in EM fungal biomass, whereas 
the addition of N fertilizer reduced C allocation to EM 
fungi and shifted the incorporation of ^^N tracer from 
EM fungi to tree foliage. Thus, changes in C allocation to 
EM fungi in response to increased N availability could 
directly affect soil C dynamics (Clemmensen et al. 2013) 
in addition to having cascading effects on the cycling of 
N in forest soils (Hogberg et al. 2011; Nilsson et al. 2012; 
Nasholm et al. 2013). 

In addition to '^N tracer studies, the ^^N:'*N ratio 
((5'^N) of EM sporocarps and tree foliage at natural abun- 
dance levels may provide a unique opportunity to explore 
the in situ movement of N in EM symbioses (e.g.. Hobble 
and Colpaert 2003; Hobble and Hobble 2006; Hogberg 
et al. 2011). In general, EM fungi tend to be enriched in 
'^N under N-limited conditions, whereas foliage of host 
plants is depleted in '^N (Gebauer and Dietrich 1993; 
Hogberg et al. 1996, 1999; Hobble and Colpaert 2003; 
Hobble 2006; Hobble and Hobble 2008). Differences in 
^^N signatures between EM fungi and their host appear 
to be the result of fractionation of roughly 8-10 %„ dur- 
ing the creation of transfer compounds (Macko et al. 
1986; Hobble and Colpaert 2003), with preferential reten- 
tion of the heavier ^^N isotope by EM fungal biomass 
and the transfer of the lighter '*N isotope to host plants 
(Taylor et al. 1997; Hogberg et al. 1999; Kohzu et al. 
2000). Previous studies, with a phytocentric perspective, 
have used these differences in (5'^N to estimate the frac- 
tion of N in trees that enters through mycorrhizal fungi 
(Hobble and Hobbie 2006, 2008; Yano et al. 2010; Mayor 
et al. 2012; Pena and PoUe 2014). In this study, we take a 
more mycocentric perspective and evaluate whether 
changes in the <5'^N signature of EM sporocarps in 
response to N additions may be used as a relative index 
of EM fungal sink strength for N, with higher '^N values 
indicating greater amounts of N retained in EM sporo- 
carp biomass. 

This study was undertaken to examine the short-term 
(<6 years) and long-term (>40 years) effects of N addi- 
tion on EM functioning in two oligotrophic Pinus sylves- 
tris L. forests in northern Sweden. Both sites are 
characterized by multiple N addition treatments, which 
allowed us to examine the response of EM sporocarp pro- 
duction and species richness to different N application 
rates and to compare short with long-term responses. 
Additionally, after 20 and 39 years of N additions, respec- 
tively, two of the N addition treatments in the long-term 
N addition experiment were terminated, thereby provid- 
ing a unique opportunity to examine the recovery of EM 
sporocarps after the cessation of high N loading. In 
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addition to measuring cumulative EM sporocarp biomass 
and species richness, we also determined the N% and 
d^^N in EM sporocarps and tree foliage within the differ- 
ent N treatments to better understand how N additions 
affect the functional role EM fungi play in the forest N 
cycle. We hypothesize that (i) N additions will have a 
negative effect on EM sporocarp biomass and species 
richness, (ii) changes in EM sporocarp production will 
depend on the amount and duration of N additions, and 
(iii) EM sporocarp d^^N values are correlated with the 
cumulative amount of N retained in EM sporocarp bio- 
mass and decrease with increasing N availability. We fur- 
ther hypothesize that after the termination of long-term 
N additions, the abundance, richness, and <5'"'N values of 
EM sporocarps as well as foliar d^^N will become more 
similar to the control than in ongoing N addition treat- 
ments, indicating a recovery of the functional role EM 
fungi play in the forest N cycle after the cessation of high 
N loading. 

Materials and Methods 
study sites 

Two boreal Pinus sylvestris L. forests in northern Sweden 
were selected for the study. The first experimental forest 
site was located near Norrliden (64°21'N, 19°45'E) and 
was planted in 1953 after prescribed burning in 1952 
(Tamm et al. 1999). The soil is a glacial till with sand as 
the dominant fraction. The N treatments started in 1971 



and consisted of three annual addition rates of NH4NO3: 
35 (Nl), 70 (N2T), and 110 kg-N-ha"^-yr"' (N3T), 
respectively, and a control (NO) receiving ~2 kg-N-ha^^ 
yr^' by natural deposition {n = three 30 x 30 m plots 
per N treatment). The amount of N added to these plots 
has varied slightly over the years and is summarized in 
Hogberg et al. (2011). In 1991, after 20 years of N addi- 
tions, the highest N treatment (N3T) was terminated, and 
in 2009, the N2T treatment was also terminated. More 
details of soil properties among the treatments are given 
in Table 1. 

The other site was located in a homogenous 70-years- 
old P. sylvestris stand 50 km northwest of Umea at the 
Rosinedalsheden experimental forest (64°10'N, I9°45'E). 
The soil is fine sand, and the profile formed is that of a 
weakly developed podsol. At this site, three 15 ha plots 
were selected based on similar stand structure for a large- 
scale eddy covariance study examining the effects of N 
additions on net ecosystem CO2 exchange with the atmo- 
sphere. Although these N treatments were unreplicated 
due to logistical and financial constraints, they provide 
insights into ecosystems processes at a spatial scale which 
would be impossible to achieve with smaller scale, more 
easily replicated manipulations (Carpenter 1996; Sullivan 
1997; Osmond et al. 2004). Starting in 2006, each plot 
has been subjected to different N treatments: a low and 
high N addition plots that receive annual additions of 
NH4NO3 of 20 and 100 kg-N-ha^'-yr"', respectively, and 
an unfertilized control plot that receives ~2 kg-N-ha^'- 
yr^' by natural deposition. The amount of N added to 



Table 1. The cumulative amount of N added, percent N (%), 5^^f4 (°„J and C:N ratio of tine organic mor-layer, and foliar N (%) and S^^H in dif- 
ferent N addition treatments in two Pinus sylvestris boreal forests in nortlnern Sweden. At Rosinedalsheden, N addition treatments started in 2006 
and consisted of two treatments: a low (20 kg-N-ha"^-yr"') and high N additions (100 kg-N-ha^'-yr"^), and an unfertilized control. At Norrliden, 
N addition treatments began in 1971 and consistent of three annual addition rates of NH4NO3: 35 (Nl), 70 (N2T), and 110 kg-N-ha"'-yr"' (N3T), 
respectively, and a control (NO). Significant differences (P < 0.05) among N treatments are indicated by different letters. 



Site 


Nitrogen treatment 






Rosinedalsheden 


Control^ 




Low 


High 


Cumulative N added (kg-N-ha^^) 
Nitrogen (%) in the mor-layer 
■^^^N (%„) in the mor-layer 
C:N ratio in the mor-layer 
Foliar N (%) 
Foliar S'^'N (%J 


0.97 ± 0.01 
-0.99 ± 0.34 
44.83 ± 2.20 

1.13 ± 0.03 
-5.16 ± 0.90 




120 

1.17 ± 0.05 
-1.06 ± 0.34 
39.10 ± 0.54 

1.23 ± 0.07 
-4.73 ± 0.46 


600 

1.20 ± 0.02 
-0.95 ± 0.40 
38.38 ± 0.45 

2.20 ± 0.24 
-1.94 ± 0.22 


Norrliden 


NO^ 


Nl 


N2T 


N3T 


Cumulative N added (kg-N-ha"') 
Nitrogen (%) in the mor-layer 
^^^N {%„) in the mor-layer 
C:N ratio in the mor-layer 
Foliar N (%) 
Foliar 5'=N (%„) 


0.96 ± 0.03" 
-0.01 ± 0.25 
45.31 ± 1.33" 

1.21 ± 0.01" 
-3.13 ± 0.13" 


1320 

1.45 ± 0.09*= 
-0.83 ± 0.17 
33.54 ± 2.42'' 

1.42 ± 0.03*= 
-0.85 ± 0.25*= 


2460 

1.58 ± 0.07'' 
-0.16 ± 0.13 
28.64 ± 0.51'' 
1.54 ± 0.02" 
0.44 ±0.14" 


2070 

1.41 ± 0.11'' 
-0.07 ± 0.17 
30.61 ± 0.22'' 

1.38 ± 0.03'' 
-2.06 ± 0.40" 



'Control treatments at both sites received only ambient deposition of N (~2 kg-N-ha '-yr '). 
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the high N plot is typical of many scientific N addition 
experiments (see Hasselquist et al. 2012) and is used here 
to indicate the sensitivity of the system to N additions. In 
contrast, the dose in the low N plot was designed to 
mimic potential future N deposition rates in southern 
Scandinavia according to model projections (Lamarque 
et al. 2005). More details of soil properties among the 
treatments are given in Table 1. 

Ectomycorrhizal sporocarp and tree foliage 
measurements 

Within each of the N treatment plots (30 x 30 m) at 
Norrliden, a 10 x 10 m plot was installed for the collec- 
tion of EM sporocarps. At Rosinedalsheden, we installed 
three 10 x 10 m plots within each of the three 15 ha N 
treatment plots. Thus, at both forest sites, EM sporocarps 
were collected from three 100 m^ plots per N treatment. 
Plots were sampled weekly throughout the summer, but 
most samples were collected between August and Septem- 
ber, as this period corresponds with most EM sporocarp 
emergence. All fully developed EM sporocarps were col- 
lected and brought back to the laboratory where they 
were identified. At the end of the growing season, cur- 
rent-year needles were collected from the top whorls of 
10 trees within each 100 m^ plot at both sites for nutrient 
and isotopic analyses. Additionally, five soil samples were 
collected in each plot using a 2.5-cm-diameter soil corer 
for nutrient and isotopic analyses of the organic mor- 
layer. 

Preparation and sample analyses 

Samples were dried (70°C, 48 h), weighed, and then 
ground in a ball mill (MM200, Retsch GmbH, Haan, Ger- 
many). Thus, dry weight and number of EM fungal 
sporocarps collected per species were recorded each week. 
If multiple sporocarps of the same EM species were col- 
lected in a plot on a given sampling date, then they were 
pooled prior to isotopic analyses. Dried samples of sporo- 
carp caps, soil mor-layer, and foliage were analyzed for 
'"■N, %N, and %C using an elemental analyzer (Flash EA 
2000) coupled to an isotope ratio mass spectrometer 
(Delta V Advantage, Thermo Scientific, Bremen, 
Germany). Results are expressed in the standard notation 
(5^^N) in parts per thousands (%„) relative to atmo- 
spheric N2. The standard deviation based on analysis of 
replicated samples was <0.1%„ for <5'^N. 

Statistical analyses 

The following community indices were calculated for each 
N treatment: 1) EM taxon richness (S), that is, the 



number of species identified, and 2) cumulative sporocarp 
abundance and biomass. Because many EM species were 
unequally represented among the N treatments and it is 
well known that EM species differ greatly in (Taylor 
et al. 2003; Hobble and Agerer 2010), this created a 
potential bias when calculating mean values for the differ- 
ent N treatments. We therefore calculated weighted means 
(i.e., community means) that took this into account. 
Community (5^^N mean values were determined separately 
for each plot by multiplying the relative abundance of 
each EM genera by the (5^^N values for that particular EM 
genus and summing these values for all EM genera found 
within a plot. Similarly, the cumulative amount of N 
retained in EM sporocarps was calculated for each plot by 
multiplying the nitrogen concentration for each EM gen- 
era by the cumulative biomass for that particular genus 
and summing these amounts for all EM genera found 
within a plot. 

At Norrliden, community indices and EM community 
(i'^'N values were subjected to a one-way ANOVA fol- 
lowed by post hoc Tukey's tests to determine significant 
{P < 0.05) differences among N treatments. For data that 
did not meet the assumption of homogeneity of variance, 
we used post hoc Games-Howell tests to determine sig- 
nificance. At each site, a two-way ANOVA followed by 
post hoc Tukey's tests was used to detect significant 
(P < 0.05) differences in (5^^N of EM sporocarps at the 
generic taxonomic level, N addition treatments, and the 
interaction between EM genera and N addition treat- 
ments. Because of the limited number of Boletus, Suillus, 
and Xerocomus taxa collected, we grouped these species 
into Boletales, which represents a higher taxonomic level, 
prior to statistical analysis. For EM genera collected at 
Norrliden, we also performed a one-way ANOVA fol- 
lowed by post hoc Tukey's tests to determine significant 
differences among N addition treatments. Given that the 
N treatments were unreplicated at Rosinedalsheden, we 
did not perform statistical analyses to determine differ- 
ences among treatments. Instead, we present treatment 
means (±SE) derived from the three plots to provide an 
estimate of within-plot variability as a guide to how 
robust the measured mean differences were among treat- 
ments. 

Because foliar N concentrations provide a useful index 
to asses soil N availability (Aerts and Chapin 2000; Hob- 
ble and Gough 2002; Cox et al. 2010), we used foHar %N 
data to examine how changes in N availability affect 
cumulative EM sporocarp biomass, species richness, and 
community ^^^N values. Least square regression was also 
used to look for significant relationships among cumula- 
tive EM sporocarp biomass, EM community (5^^N values, 
and the cumulative amount of N retained in EM sporo- 
carps. All statistical analyses were performed using SPSS 



3018 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



N. J. Hasselquist & P. Hogberg 



Response of EM Functioning to Nitrogen Additions 



Statistical software (SPSS 21.0 for windows, SPSS Inc., 
Chicago, IL). 

Results 

EM sporocarp production and richness 
among N treatments 

A total of 1657 EM sporocarps representing 11 genera 
and 33 species were collected in this study. Across both 
sites, cumulative EM sporocarp biomass and species rich- 
ness significantly decreased with increasing N availability, 
as measured by foliar N concentrations (Fig. 1). At 
Rosinedalsheden, sporocarp production and species rich- 
ness were similar between the control and low N treat- 
ment (Table 2), with Cortinarius species being the 
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Figure 1. Correlation between treatment means of foliar N and 
cumulative EM sporocarp biomass (A) and EM sporocarp richness (B) 
at Norrliden (filled symbols) and Rosinedalsheden (open symbols). 
Data points represent the mean ±SE (n = 3). 



dominant genus in both treatments (Fig. 2A). In contrast, 
only 12 EM sporocarps were collected in the high N treat- 
ment, and these were found in only one plot. Moreover, 
aU 12 sporocarps in the high N treatment were Laccaria 
bicolor (Maire) Orton, whereas L. bicolor was not found 
in the other two treatments (Fig. 2A). A complete list of 
EM sporocarps found among the N treatments at Rosine- 
dalsheden is shown in Supporting Information Table SI. 

At Norrliden, the ongoing addition of 35 kg-N-ha^'- 
yr^' in the Nl treatment severely reduced EM sporocarp 
abundance and richness relative to the control, NO treat- 
ment (Table 2). Three years after the termination of N 
additions in the N2T treatment, EM sporocarp abundance 
and richness was still reduced compared with the control. 
Richness and cumulative EM sporocarp biomass in the 
N3T treatment was more similar to the control compared 
with the Nl and N2T treatments. In the control, roughly 
80% of all EM sporocarps collected were Cortinarius spe- 
cies (Fig. 2B). In contrast, the two most abundant EM 
sporocarp taxa in the Nl treatment were Boletales and 
Lactarius. It is worth pointing out that Boletales in the 
Nl treatment consisted entirely of Xerocomus subtomen- 
tous (L.) Quel, whereas Boletales in the other N treat- 
ments were characterized by Suillus variegatus (Swartz ex 
Fr.) O. Kuntze. The most abundant EM sporocarp taxon 
in the N2T treatment was Laccaria, whereas the recovery 
of EM species richness in the N3T treatment was largely 
the result of the re-occurrence of Cortinarius species. A 
complete list of EM sporocarps found among the N treat- 
ments at Norrliden is shown in Supporting Information 
Table S2. 

EiVI community and foliar ^^N among N 
treatments 

Across both sites, EM community 5'^N values were nega- 
tively correlated with increasing N availability (Fig. 3A), 
yet positively correlated with increasing cumulative EM 
sporocarp biomass (Fig. 3B). Higher EM community i5'^N 
values were in turn significantly correlated with greater 
cumulative N retained in EM sporocarps (Fig. 3C). 

After 6 years of N additions at Rosinedalsheden, (5'^N 
of current-year needles differed little between the control 
and low N treatment, yet were -3%^ higher in the high N 
treatment (Table 1). EM community S^^N values were 
similar between the control and low N treatment, whereas 
for the one plot in the high N treatment in which we 
were able to calculate an EM community ^'^N value, it 
was ~5%o lower compared with the control and low N 
treatment (Table 2). Taken together, differences in S^^N 
between EM sporocarps and tree foliage (A) differed little 
between the control and low N treatment, whereas A was 
much smaller in the high N treatment. 
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Table 2. Comparison of EM sporocarp taxon richness (S), cumulative EM sporocarp abundance, biomass and amount of N retained in EM sporo- 
carps, EM community C:N ratio and S^^bi, and tine '^N enrichment between EM sporocarps and tree foliage (A) among the different N treatments 
at Norrliden and Rosinedalsheden. Values shown are average ± SE (n = 3). Significant differences (P < 0.05) among N treatments are indicated 
by different letters. 



Site 


Nitrogen treatment 






Rosinedalsheden 


Control 


Low 


High 


Taxon richness (S) 

Cumulative abundance (ind per 100 m^) 
Cumulative biomass (g per 100 m^) 
Cumulative EM sporocarp N (g N per 100 m^) 
EM community C:N ratio 
EM community ^'^N (%„) 

^ ( l^sporocarps ~ ^foliage) 


13.67 ± 0.33 
162.00 ± 38.97 
78.34 ± 16.06 

2.85 ± 0.59 
12.55 ± 0.12 

5.60 ± 0.19 
10.78 ± 0.27 


12.00 ± 2.52 
193.33 ± 50.83 
90.40 ± 12.6 

3.38 ± 0.48 
12.45 ± 0.10 

5.30 ± 0.29 
10.03 ± 0.74 


0.33 ± 0.33 
4.00 ± 4.00 
0.57 ± 0.57 
0.04 ± 0.04 
7.22 
0.46 
2.49 


Norrliden 


NO N1 


N2T 


N3T 


Taxon richness (S) 

Cumulative abundance (ind per 100 m^) 
Cumulative biomass (g per 100 m^) 
Cumulative EM sporocarp N (g N per 100 m^) 
EM community C:N ratio 
EM community ^'^N (%„) 

^ ( N5porocarps ~ ^foliage) 


7.33 ± 0.88= 2.67 ± 0.67'' 
58.67 ± 17.15 13.33 ± 3.18 
52.56 ± 11.04 25.93 ± 7.35 

1.90 ± 0.40 1.40 ±0.47 
12.91 ± 0.30= 9.15 ± 0.97'' 

7.30 ±0.17= 3.98 ± 0.70'' 
10.43 ± 0.28= 4.83 ± 0.60'' 


5.00 ± 1.16=''' 
33.33 ± 14.75 
29.19 ± 15.68 

1.37 ± 0.72 
10.09 ± 0.12'' 

3.23 ± 1.02" 

2.79 ± 1.13" 


7.00 ± 1.53=''' 
87.67 ± 39.31 
54.07 ± 12.16 

2.20 ± 0.46 
11.39 ± 0.17=''' 

4.78 ± 1.23=''' 

6.84 ± 1.63= 



dependent on the amount and duration of N additions. At 
Rosinedalsheden, the annual addition of 20 kg-N ha^'-yr^^ 
for the past 6 years affected neither the abundance, 
richness, nor the cumulative amount of N retained 

in EM sporocarps, whereas the addition of 
100 kg-N-ha^'-yr^' during the same period nearly elimi- 
nated EM sporocarps. This short-term response is in con- 
trast to the long-term response observed at Norrliden, 
where annual additions of N have been going on for the 
past 40 years. Despite a considerable decrease in the abun- 
dance and a shift in the dominant EM sporocarp taxa, 
long-term addition of N at Norrliden at rates of ca. 35 or 
70 kg-N-ha^^-yr^' did not eliminate C allocation to EM 
sporocarps. Furthermore, 20 years after the termination of 
N addition in the highest N treatment (N3T), the abun- 
dance, richness, and ^^^N values of EM sporocarps and tree 
foliage were not different between the control and N3T 
treatment, indicating a recovery of the functional role EM 
sporocarps play in the forest N and C cycles after the 
cessation of long-term N loading. 

The contrasting responses of EM sporocarp communi- 
ties to the two N addition rates at Rosinedalsheden sug- 
gest that the response of EM fungi to increasing N 
availability is more complex than the often assumed con- 
stant, negative linear relationship (Arnolds 1991; Lilleskov 
et al. 2001, 2002a; Kjoller et al. 2012). The abundance 
and species richness of EM sporocarps was similar 
between the control and low N treatment, whereas the 
addition of 100 kg-N-ha^^-yr^' nearly eliminated EM 
sporocarps in the high N treatment. This finding suggests 



At Norrliden, (5 N of current-year needles was signifi- 
cantly lower in the control and N3T treatment compared 
with the Nl treatment, which was significantly lower rela- 
tive to the N2T treatment (Table 1). EM community 
(5^^N values were lowest in Nl and N2T treatments, inter- 
mediate in the N3T treatment, and highest in the control 
(Table 2). Taken together, A was greatest in the control 
and N3T treatment, followed by the Nl treatment, and 
lowest in the N2T treatment. 

^^N of EM sporocarps 

At both sites, §^^N of EM sporocarps differed significantly 
among genera (Fig. 4), whereas differences among N 
addition treatments were only marginally significant at 
Norrliden (P = 0.08). For EM taxa collected throughout 
all the N treatments at Norrliden (i.e., Boletales, Lactarius, 
and Russula), sporocarp ti'^N values were lowest in the 
Nl and N2T treatments, whereas sporocarp 5^^N in the 
N3T treatment approaching those in the control 
(Fig. 4B). In contrast, (5'^N values of Laccaria and Paxillus 
were lowest in the N3T treatment relative to the other N 
addition treatments. 

Discussion 

Results from this study clearly show a reduction in the 
amount of C allocated to and N retained in EM sporocarps 
with increasing N availability in two oligotrophic boreal 
pine forests. However, this response was strongly 
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Figure 2. Mean (±SE) relative abundance of EM sporocarp taxa at 
Rosinedalsheden (A) and Norrliden (B). Relative abundance was 
calculated for each plot as the number of each EM sporocarp taxa 
from the total number of EM sporocarps collected (n = 3 per N 
treatment). 

that with increasing amounts of N additions in N-limited 
boreal forests, there appears to be a threshold that only 
when reached elicits a negative response in EM fungal 
communities. Previously, Wallenda and Kottke (1998) 
suggested a critical load of 35 kg-N-ha^'-yr^' as sufficient 
to cause changes in the abundance of EM sporocarps, 
which is consistent with our observation of no negative 
effect of 20 kg-N-ha^^-yr^', but a strong negative effect of 
100 kg-N-ha^' yr^^ at Rosinedalsheden. 

Despite a considerable change in EM sporocarp taxa, 
long-term additions of 35 and 70 kg-N^^-ha^'-yr^' in 
both the Nl and N2T treatments; respectively, at Norrl- 
iden did not eliminate C allocation to EM sporocarps, 
although the amount of C allocated to EM sporocarps 
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Figure 3. Correlation between treatment means of foliar N and EM 
community iS'^N (A), cumulative EM sporocarp biomass and EM 
community (B), and EM community and the cumulative 
amount of N retained in EM sporocarps (C) at Norrliden (filled 
symbols) and Rosinedalsheden (open symbols). Data points represent 
the mean ±SE (n = 3). 

was small (<100 g-m^^). Increasing N additions resulted 
in higher dominance of Lactarius and Laccaria species 
(e.g., nitrophilic species) in the Nl and N2T treatments, 
which are presumably more successful at using inorganic 
N sources, have hydrophilic short- and medium-distance 
or contact exploration types, and thus pose a lower C cost 
to their host (Hobble and Agerer 2010; LiUeskov et al. 
2011). Cortinarius species, on the other hand, were the 
dominant EM taxa in the control, yet absent from the Nl 
and N2T treatments. Cortinarius species are thought to be 
more capable of using complex organic N sources 
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Figure 4. Mean (±SE) i5^^N values of EM genera found within the 
different N addition treatments at Rosinedalsheden (A) and Norrliden 
(B). Different capitai letters at the top of each graph represent 
significant (P < 0.05) differences among EM taxa. For each EM taxa at 
Norrliden, F and P values (with df) are from a one-way ANOVA 
testing for the effect of N addition with significant differences 
indicated by different lower-case letters. 

(Lilleskov et al. 2002b), possess hydrophobic long- 
distance exploration types (Hobble and Agerer 2010), and 
consequently may require more C from their host (Lilles- 
kov et al. 2011). The notion that changes in EM fungal 
communities in response to N additions at Norrliden 
may be the result of taxon-specific differences in N use is 
consistent with a previous study at the same site that 
showed significantly higher concentrations of inorganic N 
in the Nl treatments compared with the control 
(Hogberg et al. 2006). However, these broad generaliza- 
tions about the source of N used by EM fungi are largely 
based on laboratory studies and remain untested in the 
field. 

The termination of two long-term N addition treat- 
ments at Norrliden provided a unique opportunity to 
examine the temporal recovery of EM sporocarp produc- 
tion after N saturation. Three years after the termination 
of N addition in the N2T treatment, the abundance and 
richness of EM sporocarps was still reduced compared to 



the control, indicating no immediate recovery of EM spo- 
rocarp production. In contrast, EM sporocarp abundance 
and richness in the N3T treatment was similar to the con- 
trol, suggesting that within 20 years after the cessation of 
long-term N additions, C allocation to EM sporocarp 
production was restored. This finding is in contrast to 
previous studies that have shown the production of EM 
sporocarps can be depressed for decades or become extir- 
pated after high N loading (Strengbom et al. 2001; Ar- 
nolds 2010). The recovery of EM sporocarps in the N3T 
treatment was largely the result of the re-occurrence of 
Cortinarius species, which provides further support that 
Cortinarius sporocarps may be a good indicator of N 
availability in N-limited boreal forest (Brandrud 1995). 
However, it is worth pointing out that Lactarius species 
were also abundant in the N3T treatment suggesting a 
long-lasting legacy effect of N additions on EM sporocarp 
richness. 

Results from this study provide evidence for a tight 
coupling between tree C allocation to and N retention in 
EM sporocarps and moreover highlight the potential use 
of (5'=*N in EM sporocarps as a relative index of EM fun- 
gal sink strength for N. In sites characterized by low N 
availability, trees relied heavily on EM fungi for N uptake, 
as indicated by lower foliar ^^^N values and greater differ- 
ences in (5'^N between EM sporocarps and tree foliage 
(A) (Taylor et al. 2003; Hobble and Hobble 2006; 
Hogberg et al. 2011). In exchange for improving N 
uptake, trees allocated more C to EM sporocarps, which 
resulted in greater amounts of N retained in EM sporo- 
carps, because of the high demand for N to build fungal 
biomass, and higher EM community (i'^'N values. How- 
ever, with increasing N availability, trees become less 
dependent on mycorrhizal fungi for N uptake, as indi- 
cated by a reduction in EM root tips (KjoUer et al. 2012), 
and as a result, N is readily taken up by tree roots. There 
is little fractionation of N when directly taken up by plant 
roots (Evans 2001), which in part may explain similar 
foliar and soil i5'^N values in the ongoing Nl treatment. 
Additionally, less N transferred by EM fungi would result 
in reduced N fractionation by EM fungi, which in turn 
could explain the negative relationship between EM com- 
munity (5^^N values and increasing N availability. Thus, 
changes in EM community (5'^'N values in response to N 
additions appear to be the result of less C allocated to 
and correspondingly less N retained in EM sporocarps 
and thus may provide a relative index of N retention by 
EM fungi. 

We recognize that differences in EM community ^^^N 
values among treatments could also reflect changes in 
dominant EM genera and the depth at which their myce- 
lium occurs. It is well known that ^'^N values increase 
with soil depth (Nadelhoffer and Fry 1988; Hogberg 
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1997), and more recently, Hobbie et al. (2014) demon- 
strated that Cortinarius species use soil N at deeper 
depths than Laccaria and Lactarius species. The observed 
differences in EM community (5^^N values among N addi- 
tion treatments may therefore simply reflect changes in 
dominant EM taxa and where they access N in the soil 
profile. However, sporocarp i5'^N for individual EM taxa 
also differed among N addition treatments at Norrliden, 
and moreover, these differences were positively correlated 
with changes in their abundance within the different N 
treatments. Taken together, our results suggest that 
changes in EM community (5^^N values in response to 
increasing N availability were not entirely the result of 
changes in dominant EM sporocarp taxa, but also reflect 
the demand for N to build fungal biomass (e.g., sporo- 
carps) which is strongly dependent on the amount of C 
allocated to EM fungi. 

Interestingly, 20 years after the termination of N addi- 
tions in the highest N treatment (N3T), EM community 
^'^N and A values approached those observed in the con- 
trol. Furthermore, when expressed as a percentage of the 
total amount of soil N, the cumulative amount of N 
retained in EM sporocarps was twice as great in the control 
and N3 treatment compared with the Nl and N2T treat- 
ments. Although the total amount of N retained in EM 
sporocarp is small, ranging between 0.14 and 
0.22 kg-N-ha^'-yr^', it is important to point out that EM 
sporocarps represent only a small fraction of the total EM 
fungal biomass in natural ecosystems (Sims et al. 2007), 
and therefore, a much larger portion of N retained in EM 
fungal biomass would be found in extramatrical mycelium 
in the soil (Wallander et al. 2004; Nilsson et al. 2012). 
Indeed, it has been estimated that 2.5-7.5 kg-N-ha^^ can 
be retained in EM mycelium each year, which represents a 
significant portion of the total amount of soil N retained 
in European coniferous forests (Nilsson et al. 2012). In 
addition to representing a large pool of organic N, EM 
fungal litter is thought to mineralize slowly and so could 
be important for ecosystem N retention (Langley et al. 
2006). Taken together, these findings suggest a recovery of 
EM functioning after cessation of high N loading and also 
highlight the potential role EM fungi play in N retention 
in forest soils (Hogberg et al. 2011; Nilsson et al. 2012; 
Nasholm et al. 2013). However, further studies aimed at 
determining the decomposition rates of EM fungal litter in 
the field are needed to better understand the functional 
role EM fungi play in the forest soil N cycle. 

In this study, we choose to focus on EM sporocarps 
when evaluating the effects of N additions of EM fungi 
because of their known sensitivity to N availability as well 
as the relative ease to identify and perform isotopic analy- 
ses on individual species. We recognize that EM sporo- 
carps represent only a fraction of total EM biomass (Sims 
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et al. 2007) and that changes aboveground do not neces- 
sarily correspond to changes belowground (LiUeskov et al. 
2002a,b; Avis et al. 2003; Toljander et al. 2006). Neverthe- 
less, our results are consistent with other studies that have 
shown a reduction in C allocation to mycorrhizal fungi in 
response to N additions at high dosage rates or over long 
periods of time (Brandrud 1995; Wallenda and Kottke 
1998; LiUeskov et al. 2002a,b; Avis et al. 2003; Toljander 
et al. 2006; Parrent and Vilgalys 2007; LiUeskov et al. 
2011; KjoUer et al. 2012). Additionally, the strong correla- 
tions among cumulative EM sporocarp biomass, EM 
community and the amount of N retained in EM 

sporocarps provide evidence that the production and 
(5'^N of EM sporocarps could be used as relative indica- 
tors of the functional role EM fungi play in the forest C 
and N cycles in light of increasing N availability. 

Changes in C aUocation to EM fungi in response to N 
additions could also alter aboveground plant growth 
through a process known as the mycorrhizal-fungal- 
induced progressive limitation (PNL) hypothesis (Alberton 
et al. 2007). According to this hypothesis, greater below- 
ground tree C allocation in N-limited environments is 
coupled to a higher EM fungal sink strength for N, and 
consequently, less N is transferred to host plants resulting 
in a reduction in aboveground plant growth (Nylund and 
WaUander 1989; Colpaert et al. 1992, 1996; Alberton et al. 
2007; Correa et al. 2008). The negative relationships 
among nitrogen availability, C aUocation to and N reten- 
tion in EM sporocarp biomass, and foliar N concentrations 
found in this study support this hypothesis and could help 
explain reduced aboveground tree growth in the control 
relative to the Nl and N2T treatments as previously 
described in Hogberg et al. (2006). Interestingly, the recov- 
ery of C allocation to EM sporocarps in the N3T treatment 
may provide an explanation for the observed reduction in 
aboveground tree growth in the N3T treatment after the 
termination of N additions (Hogberg et al. 2006). 

Conclusions 

This study clearly demonstrates that the response of EM 
sporocarps to increasing N availability is complex and is 
strongly dependent on the amount and duration of N 
additions. At Rosinedalsheden, the addition of 
100 kg-N-ha^'-yr^^ for 6 years essentially eliminated EM 
sporocarps, whereas the addition of 20 kg-N-ha^'-yr^^ 
did not affect sporocarp production and if anything 
slightly increased the amount of N retained in EM sporo- 
carps. Although long-term additions of N at Norrliden 
lead to a change in the dominant EM sporocarp taxa, it 
did not eliminate EM sporocarp production. Instead, 
reduced tree C allocation to EM sporocarps in response 
to long-term N additions was coupled to lower EM 
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community <5'^N values, suggesting a decrease in the 
functional role of EM fungi as N sinks. However, 20 years 
after the termination of N additions in the most heavily 
fertilized treatment, EM sporocarp abundance, richness, 
community i5^^N, and A values were no longer different 
than the control, suggesting a recovery of the functional 
role EM fungi play in the forest N cycle. Given the 
important functional role EM fungi play in forest soil C 
and N dynamics (Hogberg et al. 2011; Nilsson et al. 2012; 
Clemmensen et al. 2013; Nasholm et al. 2013), a better 
understanding of the dosage-dependent effects and recov- 
ery of EM functioning after N saturation should be a 
strong future research priority. 
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